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Many proposed and realized spintronic devices involve spin injection and accumulation at an 
interface between a ferromagnet and a non-magnetic material. We examine the electric field, voltage 
profile, charge distribution, spin fluxes, and spin accumulation at such an interface. We include the 
effects of both screening and spin scattering. We also include both the spin-dependent chemical 
potentials and the effective magnetic field H* that is zero in equilibrium. For a Co/Cu interface, 
we find that the spin accumulation in the copper is an order of magnitude larger when both chemical 
potential and effective magnetic field are included. We also show that screening contributes to the 
spin accumulation in the ferromagnet; this contribution can be significant. 

PACS numbers: 75.70.Cn, 72.25.-b, 75.47.-m, 75.76.+] 



I. INTRODUCTION 

Although electronic current has been studied since the 
early 19*^ Century, spin current has been studied only 
much more recently. In particular, spin transport across 
interfaces between metals and ferromagnets has been an 
important topic sin ce the discovery of giant magnetore- 
sistance (GMR),'32]the principle behind the predominant 
method of reading stored data. The magnetic read-head 
of a hard drive contains a thin non-magnetic layer sand- 
wiched between two ferromagnetic layers. 

A theory for spin current and electrical potential at 
a metal/ ferromagnet interface is given by Johnson and 
Silsbee^ (JS); an appendix of that work is devoted to 
electric currents crossing such interfaces, and it consid- 
ers the effect of the interfaces on spin fluxes and on elec- 
trical voltage. Detailed theories for electrical currents 
crossing metal/ferromagnet multilayers (that is, series of 
interfaces) are given by Valet and FerlP (VF) , which in- 
cludes solutions for the electric field and spin fluxes, and 
by Hershfield and Zhao^ (HZ). However, none of these 
theories considers semiconductors, and each makes a dif- 
ferent assumption, not made by the present work, about 
some part of the magnetoelectrochemical potential (first 
defined by JS and discussed in detail below). JS neglects 
the chemical potentials fj,^ and /z^, HZ neglects the effec- 
tive magnetic field^ H* (discussed below), and VF takes 
the chemical potential to be spin-independent. 

The present work revisits the problem of spin trans- 
port across the interface between a non-magnetic mate- 
rial (NM) and a ferromagnet (FM), and calculates the 
electric field, voltage, charge density, spin fluxes, and 
spin accumulation. The results also apply to FM/FM 
and NM/NM interfaces. We show that inclusion of both 
H* and fi^ and fi^ are necessary to predict the spin ac- 
cumulation near the interface. For copper, neglecting 
either contribution decreases the spin accumulation by 
about a factor of ten. Further, this work includes the 
surface screening mode (called the charge mode by HZ), 
neglected by JS and VF, and ultimately neglected by HZ, 
which for large screening lengths (semiconductors) plays 



an important role in determining the spin current and the 
spin accumulation. Including the screening mode permits 
the electric field and potential to be continuous across the 
interface. Previous works allow the field and potential to 
be discontinuous. Reference [71 which extends VF by cal- 
culating the spin accumulation when the non-magnetic 
material is semiconducting, also neglects screening.El 

Section [H] briefly discusses the equations that govern 
spin-dependent transport in solids. Section [lll| finds the 
deviations from equilibrium due to the screening mode, 
the spin-diffusion mode, and a bulk response associated 



with the applied electric current. Section IV discusses the 
bulk and boundary conditions at an isolated interface. 
Section [V| compares the assumptions of the present work 



to those of previous theories. Section VI gives the electric 



field, voltage, charge density, spin fluxes, and spin accu- 
mulation near a Co/Cu interface. Section VII provides 
a brief summary and conclusion. Appendix |A| shows de- 
tailed calculations for the spin-diffusion mode, the results 



of which are given in Sec. HI and App . [B] explicitly gives 
the boundary conditions discussed in Sec. |IV| 



II. TRANSPORT EQUATIONS 

We use superscripts I and II or NM and FM to denote 
adjacent materials. When developing bulk equations that 
apply separately within each material, we omit the su- 
perscript, and reintroduce it when discussing materials 
in contact (or when discussing properties specific to a 
FM or a NM). 



A. Fundamental Relations 

Within each material, the nu mbe r and current densi- 
ties n^^i and j^^^^ are related bjJ^Hni 



dtn-f + dij^^ = S, dtn^ + dij^, = -S. 



(1) 



Here S is the rate at which down-spins flip to up-spins. 
We consider the total electric current density J = — ejtot 



2 



to be a known uniform constant, and continuous across 
an interface. For current along x across an isolated in- 
terface (in the yz-plane) between materials I and II, we 

have jtot = jf j + jf j - df^ + df^ ■ 

We take M, the direction of the magnetization M, to 
be fixed. Since the electron g-factor is negative, for ma- 
jority carriers defined to have up-spins, then M is aligned 
with the down-spins. 

We take /2 to be the maqnetoelectroche mical potential, 
defined for up- and down-spin electrons a^^Hni 



fi^ = fi^ - e(l) -\ — fj.oH ■ M, 



= Mi 
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fioH* ■ M, 



(2) 
(3) 



Here a > (with units of a density of states per second) is 
related to a characteristic spin-flip time (or, equivalently, 
to a characteristic spin- flip length). 

We are interested in steady-state solutions, so that 
dtn^ = = 9tni . Taking the gradient of Eq. (7| and 
employing Eqs. ([l]) and ([s]) then gives two coupled differ- 
ential equations for /l^ and /l^. 



e 



2 Mr = 



(9) 
(10) 



On applying appropriate boundary conditions, Eqs. ^ 
and (10) give p,^ and /i^. 



where /i^ and /i^ are the respective chemical potentials 
of up- and down- spin carriers, e > is the magnitude 
of the electron charge, (j) is the electrical potential, g is 
the dimensionless g-factor (with jf/j 2 for electrons), 
/zb is the Bohr magneton (with units of J/T), and /io 
is the permeability of free spac^I^ (with units of N/A^). 
In the simplest case, H* is the difference of the exter- 
nal field Hq and the uniform exchange field Hint (with 
Hint li M). More generally, in addition to Hq we must 
include the magnetic dipole field Hdip, the crystalline 
anisotropy field Ham a-nd the non-uniform exchange field 



Hex (proportional to V Af) 



H* — Hq + Hdip + Ha 



H, 



Hi, 



(4) 



We have if* = in equilibrium.^^ 

By irreversible thermodynamics (see, for example, the 
general treatments in Refs. [TH [151 and [TBI or the spin- 
related treatments of Refs. jHI and [TUI, the non-negativity 
of the rate of entropy production implies that the fluxes 
can be written in terms of thermodynamic forces, i.e., 
gradients of intensive thermodynamic quantities. Thus, 



(5) 
(6) 



B. Linearized Relations 

We are interested not only in ft^ and fii, but also in 
and - in particular, the difference of their devia- 
tions from equilibrium Sn^ — Sn^, i.e., the spin accumu- 
lation (which is proportional to the "out-of-equilibrium 
magnetization" or "nonequilibrium magnetization" dis- 
cussed by VF and HZ). Near equilibrium, we can lin- 
earize the deviations (denoted by S) from equilibrium of 
the chemical and magnetic contributions to the magne- 
toelectrochemical potentials: the chemical potential de- 
viations can be written as 



(11) 



and the deviation in the effective magnetic field at fixed 
Hq can be written as 



5H* ■ M 



X 



M 



X 



(Jn^-Jn;), (12) 



where x is the magnetic susceptibility for an isotropic 
material (defined by Xij = X^ij)- Thus Eqs. ^ and ^ 
give 



5n 



(T,i) 



6(50 ± 



{5n^ - Sn^) 
2N^ 



(13) 



where cr^ and cr^ are the respective electrical conductiv- 
ities of electrons of up- and down- spin, and the coeffi- 
cients = L^\^ by the Onsager principle. We have im- 
plicitly neglected temperature gradients, which can also 
contribute to spin fluxes .'^^''^ Neglecting the off-diagonal 
coefficients L^^ — L^i, we have 



Jr. 



(7) 



T he n on-negativity of the rate of entropy production 



-a (/i^ - /X;) 



(8) 



where we define 



t 



X 



miJ-BiJ'O 



(14) 



each of which has units of a density of states. 

There are thus three unknowns {5n^, 5n^, and Scf)). 
Eqs. ([9| and ( 10 ) give two coupled differential equations, 
and Gauss's law provides a third: 



eoe 



{5n^ + 5n^) . 



(15) 



For the bulk response and each of the surface mode, we 
must find (5n^, (5n^, and 5(j)- 
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III. STATIC BULK RESPONSE AND SURFACE 
MODES 

We now study the static bulk response and surface 
modes of the system. For brevity we write surface so- 
lutions to have the form e~^/^ where £ is some length, 
although for the material on the left side of the interface 
one should use (because the deviations must decay as 
X — > —oo). In general, each surface solution has the form 
^±ixLnt-x)/i T^i^Qj-Q 3;.^^ is i]^Q position of the interface, but 
we take the interface to be at x-mt — 0. 

The electric field E and voltage are continuous every- 
where. (We call these "Maxwell conditions.") To ensure 
this, we include the surface screening mode. JS, VF, and 
HZ neglect screening and do not satisfy these conditions. 

We first discuss the bulk response associated with the 
electric current, which has a simpler structure than the 
surface modes associated with screening and with spin- 
diffusion. 



A. Bulk Respone (dc) 

We consider a system with a uniform constant electric 
current. The (bulk) response associated with this cur- 
rent, which can be thought of as a "dc mode" (dc), is 
characterized by a constant uniform electric field (which 
in principle differs for each material) . We define this field 
as 



Ode 



where Eq^^ is a constant determined by applying bound- 
ary conditions. The potential associated with this mode 
is 



Ode 



Eo,.x + Vr 



Ode J 



(17) 



where Vq^^ is another constant (with units of V) deter- 
mined by applying boundary conditions. By Gauss's Law 
there is no overall (bulk or surface) charge associated 
with this mode, as expected. Further, 

%dc = %dc = ~eS(l)dc = eEo^^x - eVo^^. (18) 
Equation Q gives 



JTd 



^T^Odc 



(19) 



Because does not necessarily equal (e.g., as for 
ferromagnets), there may be a non-zero spin current as- 
sociated with the dc mode. 



B. Screening Mode (Q) 



One solution to Eqs. (|9]), (10) and (15) has 5^^ = = 
(5/2^ so that = = ii- This mode is therefore entirely 



static (neither spin current nor charge current), corre- 
sponding to electric screening and characterized only by 
charge and potential gradients. We therefore designate it 
the "screening mode," and use the subscript Q to denote 
its properties.l^il 

Note that for metals the screening mode is not well- 
described by the present type of theoix but is instead 
associated with Friedel oscillations.'22H2l! xhe following 
treatment of screening is more appropriate for doped 
semiconductors. 



By Eq. (13 1, setting ^A'Tq ~ ^^^iQ =0 — = relates 
the up- and down-spin concentrations, 

N^{N^ + N^) 5n^^ = {N^ + ) 5n^^ , (20) 

and setting (5/2^^ -I- Sfl^^^ = + = Osu gives 

1 Sn:,^' 



2e V 



Define 



(21) 



(22) 



with units of a density of states. Substitution of Eq. ( 20 ) 



into Eq. (21) then yields 



(23) 



Substitution of Eqs. (20 1 and (23) into Gauss's Law, 
Eq. ( 15 ) then gives 



2e 



(16) N^dj6n. = — [N^ (iV; + N^) + {N^ + N^)] 5n. 



With the definitions 

^ ~ 2e2 Nj ' 



Eq. ( 24 ) can be written as 



(24) 

(25) 
(26) 

(27) 



For X 00 and e 1, Eq. ([25| gives £q = £o/[e^(^t + 
N^)], which agrees with Ref. P, 

We now define the quantity Vqq such that 



6n 



c 



which satisfies Eq. ( 27 1 . Then Eq. ( 20 1 gives 



Sn 



2e 



(7Vt + ^x)^o«e--/'«, 



Ar2 



(l2 "^OqC , 



(28) 



(29) 



(30) 
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and Eq. (21 ) gives 



given by 



(31) 



The screening mode can lead to a nonzero spin accu- 
mulation, defined by 



(32) 



Equations (28 1 and (29) give 



An, -^^ ^ 



2^{N^-N^)eVo^e--^'^, (33) 



NsN^ - 2Nj 



\ NsN^ ~ 2Nj 

[ NsN^ - 2Nj 



(37) 
(38) 

(39) 



which is nonzero in a ferromagnet, where ^ N^. 



C. Spin Mode (5) 



The second solution to Eqs. (|9|, (10 1 and (15) is more 
complicated than the screening mode. It is characterized 
by a nonzero spin current = j-t- — ^ 0. We therefore 
designate it the "spin mode," and use the subscript 5* to 
denote it P Following Ref. IHIwe also use Q (for charge) to 
denote the screening mode. (The reader is thus warned 
that S refers to spin, not to screening.) 

We now give the solution for the characteristic length, 
the spin concentrations, the electrical potential, and the 
spin accumulation associated with this mode. The details 
of the analysis are given in Appendix \K\ 

Define the up- and down-spin associated lengths 
and , which satisfy 



/2 ^ 



(34) 



The decay length associated with the spin mode, vari- 
ously called the "spin- 
is then given bjElEll 



ously called the "spin -flip" or "spin-diffusion" length, 4f , 



1 

w 



(35) 



For a non-magnetic material, the dimensionless coeffi- 
cient ^ — > 0. 

The spin mode leads to a nonzero spin accumulation; 
Equations (38) and (39) give 



, (40) 



I NsN^ - 2Nj 

so that Arias is nonzero for both ferromagnets and non- 



magnetic materials. For the latter, Eq. (40) simplifies 
to 



(NM) _ iVf^x'^^Os „-x/t,i 



(41) 



The spin-carrier currents associated with the spin 
mode are given by 



Vo 



(42) 



The total electric current —edjtot — ^^i^j^ + Sji) ~ 
for the spin mode, but there is a nonzero spin current 
Sja = Sj^ - 6ji, given by 



6ja 



(43) 



We also define 



D. Description Near Interface 



C: 



4s 



(36) 



where Ns has units of a density of states and C is di- 
mensionless. With V05 a constant to be determined by 
boundary conditions, the deviations in the electrical po- 
tential and up- and down- spin concentrations are then 



A full description of the region near an interface in- 
volves the combination of both surface modes (S and Q) 
derived above, and the bulk constant current (dc) mode. 
For the potential, electric field, charge density near an 
interface located at x = a;int, from Eqs. (37), (A14), (30), 



(31), ilGl, and (17) we have, with four unknowns per 



material {Eq^^, Vq^^, Vqq, and Vo,) to be determined by 
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boundary conditions, 

~EoAx-Xint) + Vo,^, (44) 

(45) 



Sp. 



-eoe 



Os ± 



(46) 



The top (bottom) sign corresponds to the material on 
the left (right) of the interface. 

The contributions to the total electric current from the 



surface modes is zero, as expected, so that Eq. ( 19 ) gives 
the electric current to be everywhere given by 



J = -ejtot = - (cr^ + CT;) i^Odc- 



(47) 



The spin mode does contribute to the nonconserved spin- 
up, spin-down, and total spin currents, which, combining 
Eqs. (42) and (19), are given by 



Jt = ■ 
Jl = - 

ja = 



Vr 



Os ±(x-a;i„t)/4t 



Ode 



^Os ^±(x-a;int)/4f 



+ / e^sf 



(48) 



(49) 



(50) 



There is no contribution from Vqq because there are no 
carrier currents associated with the charge mode. 



For the spin accumulation, Eqs. (41), (40), and (33) 
yield 



+ 2e-^(iVt - A^;)V^o«e±(---'-)/^«. 



For a non-magnetic material this simplifies to 



^^(NM) ^ ^T^xe^Os ^±{x-x,^t)/e. 



(51) 



(52) 



IV. BOUNDARY AND BULK CONDITIONS 

For an isolated interface at Xi^t = (see Fig. [T]) be- 
tween materials I (at x < 0) and II (at a; > 0), in general 
there are eight unknowns {Eq^^, ^0^^' ^Oq, and Vqs for 
each of materials I and II). There are eight conditions?23 

(i-ii): the potential and field E are continuous across 
the interface ~ Maxwell conditions; 



z Material I (FM) 



Material II (NM) 







FIG. 1. An isolated interface between a ferromagnet (dark 
gray, at a; < 0) and a non-magnetic material (light gray, at 
a; > 0). This work considers an electric current density Jx, 
and magnetization of the FM along ±£. 



(iii) : the electric current — e(j^-l-j^) is continuous across 

the interface - charge conservation; 

(iv) : the spin current is assumed continuous across the 

interface (although we take both up- and down- 
spin currents to be continuous, this is only a single 
condition since condition (iii) constrains their sum) 
- assumption of no surface spin-scattering; 

(v-vi): the up- and down-spin currents across the in- 
terface are directly proportional to the discontinu- 
ity in up- and down-spin m agnetoel ectrochemical 
potential across the interface^EHnESl _ irreversible 
thermodynamics; 

(vii) : the total electric current —e{j^ + Ji) is a known 

constant; and 

(viii) : there is an arbitrary constant voltage (which we 
define by setting the voltage V^^^ = 0). 

Conditions (i)-(vi) are boundary conditions and (vii) 
and (viii) are bulk conditions. They are explicitly cal- 
culated in Appendix [B] 

For a multilayer (a series of k interfaces between k + 1 
materials) , each additional interface adds another of each 
of the boundary conditions (i)-(vi), so that in general 
there are 6A; -|- 2 conditions. 



V. COMPARISON TO PREVIOUS THEORIES 

As noted above, the theories of Johnson and Silsbee 
(JS), Valet and Fert (VF) and Hershfield and Zhao (HZ) 
neglect the screening mode, and therefore cannot have 
field and potential continuity at the interface. Further, 
JS neglects the chemical potentials and fi^ and HZ 
neglects the internal magnetic field H* . The discrepancy 
between predicted spin accumulation, found below, par- 
ticularly in a non-magnetic material, demonstrates that 
inclusion of all parts of the magnetoelectrochemical po- 
tential is essential for calculating the spin accumulation 
in a non-magnetic material, even to within an order of 
magnitude. For comparison of Eq. (41) to the spin ac- 
cumulation predicted for these other works W = HZ and 
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W = HZ, we define the dimensionless factor C„ as 

= CwA<^). (53) 



Note that — I for the present work. We show be- 
low that if one of or is near unity (and therefore 
agrees with the present work), then the other diverges, 
so that at least one of the assumptions gives results that 
significantly disagree with the present work. 



A. Neglecting H* and the Screening Mode 

Neglecting the last term (proportional to H* ■ M) in 
Eq. (13), as in HZ,l^is equivalent to taking °o (and 
therefore — )■ oo) in the present results. Under this 



assumption, 2N^ and iV| -J> N^{N^ + N^)- Equa- 



tions (37l-(39) then simplify to 



2 (Ns -N^- N^) 



Vb„e" 



-2N^ + Nsil + C) 
2 {Ns -N^- N^) 

2N^ + Nsi-l + C) 
2 {Ns ~N^- N^) 



(54) 
(55) 
(56) 



HZ neglect the screening mode, so the spin-diffusion 
mode is the only surface mode, and it gives a spin ac- 
cumulation of 



^AN^N^ + Ns [N^ + N^ + C{N^- N^)] 
2 {Ns - iVt - N^) 



(57) 



Direct comparison can be made to the results of the 
present work in the non-magnetic material. With 
defined by Eq. (53), we have 



Chz = (l + iVt/^x)• 



(58) 



Using Tables|I]and|H) we find Chz ~ 0.0986 for Cu. (Cu is 
a diamagnet, therefore it has N^ < 0; for a paramagnet, 
where N^ > 0, the underestimation of spin accumula- 
tion for the HZ assumptions is less striking, although it 
remains significant.) For the ferromagnet, the spin ac- 
cumulation due to the screening mode is neglected, and 
the spin accumulation due to the spin mode agrees with 
the present work to within the precision of the present 
calculations. Hence, the assumptions made by HZ seem 
appropriate for ferromagnets but not for nonmagnetic 
materials. 



B. Neglecting fi^, and the Screening Mode 

JS neglects the chemical potentials fi^ and fi^ in 
Eq. (13), which is equivalent to taking N-^^i — >■ oo in 



TABLE I. Bulk and interfacial properties of cobalt and cop- 
per, and well-known constants. Here, A is the area of the 
interface, and R is the spin-dependent interface resistance. 
^Value is for the (100) orientation. *The susceptibility of 
Cobalt is field-dependent, with 70 < x'^° < 250 (see Table 2.2 
of Ref. I29[) : we take an intermediate value. 



Quantity 


Value 


Units 


Ref 


O 


2.47 X lO'^ 




1301 


''I 


0.913 X 10^ 




1301 




8.35 X 10^ 




1301 


£ ( 
^st 


59 X 10^^ 


m 


[30] 




450 X 10~^ 


m 


[30] 




5.10 X 10*" 




1301 




19.7 X 10'*'' 




1301 




3.89 X 10*'' 




1301 


^^Cu/Co 


0.31 X 10"*'^ 




mi 


^^Cu/Co 


2.31 X 10"*^ 




mi 




100 




nut 


X^" 


-0.932 X 10"^ 






Ms 


9.27 X 10"^* 


J-T-i 




MO 


4tt X lO"'^ 


N - A-2 




eo 


8.85 X 10~*^ 


A-s-V-^-m" 


1 


e 


1.6 X 10"*^ 


C 






2 







the present work. It also neglects the screening mode. 
Various properties of the spin mode are now calculated 
under these assumptions. 



Equation (37 1 gives 



c 



Further, Eqs. ^ and ^ give 

so that the spin accumulation is given by 



(59) 

(60) 
(61) 

(62) 



Direct comparison can be made to the results of the 
present work in the non-magnetic material. With C.^ 
defined by Eq. ([53|), we have 



C.,s =(l + ^xM)- 



(63) 



Using Tables|T]and|TTJ we find C,,s ~ -0.109 for Cu. Thus, 
the JS assumptions seem inappropriate for determining 
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TABLE II. Bulk and interfacial properties of cobalt and cop- 
per, calculated from the results o f th e present work (and Ta- 
ble [l| . Here ol is found from Eq. ( A4 1 . 



Quantity 


Value 


Units 


„Cu _ „Cu 1 „Cu 


16.7 X 10'' 


n-^-m 




3.23 X 10^^ 


!^"^-m 




0.433 X 10^^ 


n-^-m 


TV 


4.63 X 10^^ 


J -m 


TyCu 


—4.32 X 10 


J -m 




9.93 X 10 


J -m 


Ns 


1.71 X 10 


J -m 


ArCo 


9.26 X 10 


J -m 


AfCu 


— O.OOl X 10 


T-l ™ 

J -m 


)\rCo2 


1.15 X 10 


T-2 „ 

J -m 


2 


-3.31 X lO''^ 


J-2-m 




74.8 X 10" 


J-i-m 




8.05 X 10" 


J-i-m 


ns 


114 X 10"^ 


m 




69.0 X 10"^ 


m 


«Cu fCu 


636 X 10"^ 


m 


/iCo 

*Q 


0.0373 X 10"^ 


m 


/jCu 


0.0667 X 10"^ 


m 




-0.460 











^Co 


0.524 




^Cu 








the spin accumulation in non-magnetic materials, partic- 
ularly those that are diamagnetic. 



Note that Eqs. (58 1 and (63) preclude simultaneously 
having Cjs ~ 1 and Chz ~ 1- 



VI. CO/CU INTERFACE 

For an isolated interface (as in Fig. [T]) , Appendix |b] 
uses each of the above conditions to find an explicit equa- 
tion for the eight unknowns and writes the unknowns in 
terms of dimensionless variables. We now present numer- 
ical results for the spin fluxes (see Fig. [2]), voltage, elec- 
tric field, charge density, and spin accumulation, for a 
cobalt/copper interface, with material parameters given 
by Tables |T| and |n| 

Figures [3^-[3j; show that, outside of a screening length 
£q of the interface, the electrical potential, field and 
charge nearly coincide for the present work and HZ, with 
JS showing discrepancies near the interface in the ferro- 
magnet (a; < 0). However, the present work significantly 




-500 



500 



X (nm) 



FIG. 2. The spin-up and spin-down carrier fluxes 5j-f and 5ji 
and the total flux 5jtot and spin flux Sja, given by Eqs. (47 1 



(501, with the (uniform) total flux normalized to unity, near 
an interface between cobalt {x < 0) and copper {x > 0) vs. x 
(in nm). The grey dotted line marks 0.5, that is, half of the 
total current. 



differs from JS and HZ within a screening length of the 
interface, as seen in Figs. [3ji-[3f. Figures [Sjl and [3^ show, 
for the present work, the continuity of the electrical po- 
tential and field at the interface. They also show, for HZ 
and JS, the discontinuities in the potential and field (due 
to scale, these field discontinuities are more obvious in 
Fig. [sja than in Fig. [sjj). Figure [sj" shows, for the present 
work, the charge density due to screening. For physical 
consistency, E and must be continuous at the interface, 
so that HZ and JS must have both an infinitesimally thin 
charge layer and an infinitesimally thin dipole layer at the 
interface. 

We conclude that outside of the charge screening 
length £q (which is very short for metals), the present 
work and HZ are equally valid for calculating electrical 
potential, field, and charge, but JS differs significantly. 

Figure [4] shows the spin accumulation for the present 
work, HZ, and JS. In the non- magnetic material {x > 
0), as shown analytically in Eqs. (58) and (63), Fig. [4^ 



shows that both HZ and JS differ from the present work 
by an order of magnitude, with JS having the opposite 
sign. Fig. |4jD shows that the spin accumulation in the 
ferromagnet {x < 0) differs for the present work and HZ; 
outside of this length. Fig |4^ shows that they coincide. 
However, the spin accumulation for JS is six orders of 
magnitude larger (and not shown). This is because JS, by 
assuming that dfj.^,i/dn^^i = 0, effectively takes N^,i — 
oo so that iV^ .^ :s> N^, whereas Tables |l] and [IT] show that 
the opposite is true for cobalt. 



VII. SUMMARY & CONCLUSION 

Using irreversible thermodynamics, we predict the spin 
accumulation at an interface between two materials when 
electric current is driven across the interface. Although 
we have numerically studied a FM/NM interface, the the- 
ory also applies to FM/FM and NM/NM interfaces. 

We find that both the chemical potentials and the effec- 
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FIG. 3. The dimensionless electrical potential, field, and 



charge density (given by Eqs. (44 1- ( 46 1 ) in arbitrary units 
near an interface between cobalt < 0) and copper (x > 0) 
vs. X (in nm). Within 1 (a,b,c), and within 1 nm (d,e,f), 
of the interface. Solid line - present work; dashed line - JS; 
dotted line - HZ. HZ coincides closely with the present work 
except within a screening length of the interface. JS gives 
somewhat different results in the cobalt within a spin-diffusion 
length of the interface. 



tive magnetic field must be included to predict the spin 
accumulation in a non-magnetic material - in fact, for 
Cu the spin accumulation changes by an order of mag- 
nitude on neglect of either contribution. However, for 
ferromagnets neglecting the effective magnetic field may 
be appropriate - numerically the results are essentially 
unchanged for Co near a Co/Cu interface. 

By including the screening surface mode neglected in 
previous works, we find an additional term in the spin ac- 
cumulation for ferromagnets. For Co near a Co/Cu inter- 
face, this term decreases the spin accumulation by ^ 10% 
within a charge-screening length of the interface. Al- 
though this length is on the order of 1-10 A for metals (a 



length scale negligible in the present macroscopic theory) , 
for ferromagnetic semiconductors this length scale should 
be much larger. Note that spin injection from a fer- 
romagnetic semiconductor into a non-magnetic material 
has been observed by Refs. J3 and 34, To test this spin 
accumulation due to screening, one may apply a small 
current to an interface between, say, Ga(Mn)As and Cu. 
Using the magneto-optical Kerr effect (MOKE) , one may 
measure the magnetization (and spin-polarization) at the 
surface. We expect there to be two nonequilibrium mag- 
netization contributions near the surface, one that decays 
over the spin-diffusion length igf and one that decays over 
the screening length £q associated with screening. The 
latter effect should be more prominent in ferromagnetic 
semiconductors. 
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FIG. 4. The dimensionless spin accumulation, given by 
Eqs. (511 and (521, in arbitrary units near an interface be- 
tween cobalt (a; < 0) and copper (x > 0) vs. x (in nm). The 
spin accumulation is shown within approximately (a) 1 fim 
and (b) 1 nm of the interface. Solid line - present work; 
dashed line - JS; dotted line - HZ. In the FM (a; < 0), HZ 
nearly coincides with the present work (deviating only within 
the charge-screening length of the interface, see inset), and the 
JS- predicted spin accumulation is several orders of magnitude 
larger and not shown. In the NM, neither approximation pre- 
dicts a spin accumulation similar to the present work. 
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Appendix A: Details of the Spin Mode 

The details of the solution for the Spin Mode, whose 



results are presented in Sec. Ill C are now given 



Equations ([9])-([T0) give, with 5 denoting deviations 
from equilibrium. 



1 



Subtracting Eq. (A2) from Eq. (Al) gives 



(Al) 
(A2) 

(A3) 



where £sf is defined by Eq. (35). On neglecting SH* ■ M 



and making the identification a — (N-^/T-fi) = (N^/t^-^), 
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Eq. ([35]) agrees with Ref. HI We use Eqs. ([34]) and ([35]) 
to find a, i^^ and in terms of 4f, cr^, and cr^, since 
they are, in principle, measurable: 



(7^ 



Solving Eq. (|A3]) gives 



-x/e,t 



(A4) 



(A5) 



(A6) 



where Vq^ , with units of electric potential, is unknown, to 
be determined by boundary conditions. Since Eq. (A6) 



shows the difference in up- and down-spin magnetoelec- 
trochemical potentials to decay over the length £sf from 
an interface - this length is called the "spin-flip" or 
"spin-diffusion" length (and sometimes referred to as the 
"SDL" ) . The length igf may be measurable by employing 
the Magneto-Optical Kerr EfFect'^^ "^^ or the Inverse Spin 
Hall Effect,!^ or may be derived using GMR measure- 
ments and theory.lSH! 



/2 



Ts 
(1 



-I- cr^ 



-a;/4t 



7^ 



ft + 



(A7) 



(A8) 



Equations (|A7| and ( [AS] ) give (5/2^^ = -(f = 
— (cr^/a>)(5/i^g, which agrees with Ref. 5, Substitution of 
Eqs. ([A^ and (|A8]) into Eq.j7) gives the up- and down- 
spin carrier currents of Eq. ( 42 ) . 



We can now can write two independent relations be- 
tween (5n^ 



and 



Equations (A6) and (13) 



give the difference of the spin potentials to be 

(A9) 



and Eqs. (A7), (A8), and (13 1 give the sum of the spin 



potentials to be 

T 4- 

(AlO) 



In conjunction with Gauss's Law, Eqs. (A9) and (AlO) 



give the concentrations and electrical potential in the 
spin mode. Specifically, we use Eq. (A9) to relate Sn^^ 



solve. 



Substituting Eq. ( |All[ ) into Eq. ( |A10[ ) multiplied by 

{N^/N^){N^ + N^) gives 



+ 2e^ (N^ + N^) 



(A12) 



Substitution of Eqs. (All) and (A12) into Eq. (15) gives 



2Nj 



2 



(A13) 



The solution for S(f>s is given above as Eq. (37 1. 

Substituting Eq. ([37]) into Eqs. iAl2h and (|A11D gives 



the up- and down- spin concentrations of Eqs. (38) and 
(39). Thus, the charge distribution 5p — — e (6n^ + 5n^) 



Substitution of Eq. (|A6]) into Eqs. (JaT]) and ([A2]) yields Sps = 



associated with the spin mode is 



NsN^ - 2Nj 



iVseVo.e-^/^-', (A14) 



which is nonzero in a ferromagnet. (This result can also 
be obtained by using Eq. (37) and Gauss's Law.) Fur- 
ther, subtraction of Eq. ( 39 ) from Eq. ( 38 ) yields the spin 
accumulation of Eq. ( 40 ) . 



Appendix B: Boundary Conditions for Current 
Crossing an Isolated Interface 

Boundary conditions (i-viii) for an isolated interface 
(that is, one that is effectively an infinite distance from 
any other interface) through which an electric current 
is passed are discussed in Sec. |VI[ They are here found 
explicitly, in numerical order. 



Conditions (i-ii): From Eqs. (44) and (45), continuity 
of 54> and 5E across the interface at Xjnt = gives 

+ < + < = ^'y^ + <^ + C'. (Bi) 



^(i) 



Vr 



(I) 



^(I) 



/II) 



(11) 



^(11) 



(11) 



(B2) 



Recall that ^ = for a non-magnetic material. 

Condition (ill) : From Eq. ( 47 ) , continuity of the elec- 



to 5n<, then use Eq. (|A10| to relate dn^^ to 50s. Thus trie current across the interface gives 



Eq. (15) can be written in terms of only 5(l)s, which we 



Equation (A9) gives 



( (II) , (i: 



(11)^ ^(11) 



(B3) 



Srif 



(All) 



Condition (iv): Although the electric current is contin- 
uous everywhere, in principle at the interface there may 
be spin scattering, so that spin current is not continuous 
across the interface. However, we neglect interfacial spin 
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scattering (as is typical in this type of theory). We thus 
take 

,•(I)^n^ - .-(11)^ 



(B4) 



Using Eq. ( 48 ) , the first of these can be written as 

,(i) 




We take 



jr(o), 



(B13) 



By Eq. (B4) one may equivalently use j^^^^ = jj^^(O) and 
^ " 

-'tint 



j| '^(0). Respective substitution of Eqs. (48) and 



(B5) (|49|) into Eqs. (iBllI) and (|B12|) gives 



As discussed above, the second relation given in Eq. (B4) 
is then automatically satisfied by condition (iii), which 
constrain the sums of the up- and down-spin currents. 

Conditions (v - vi ) : The spin currents across the inter- 
face are given bySEl 



■4. ) int 1 



(B6) 
(B7) 



were, ( A)int denotes the difference between the value just 
on the right of the interface (x — 0+) and the value 
just on the left (a; 0^). Since without the electric 
field associated with the dc mode there is no steady-state 
current, the currents are proportional to the differences 
in dp, rather than p,. We now find (A (5/l)i nt for each mode 
and then substitute them into Eqs. (B6)-(B7). 

The charge mode has Sji^ = = (5/2^ , so by 



Eqs. (B6|-(B7) it does not affect the current crossing the 
boundary. At the x = x-mt = interface, Eqs. (A7) and 
(|A8l) give 
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At the interface, Eq. (17) gives 



(A^Mt<jc)int = (A(5Aii^Jint 



(I) 



(B9) 
(BIO) 



Substitution of Eqs. (|B8[)-(|B10[) into Eqs. (|B6[)-(|B7[) 

yields 



Jti: 



int 
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(B15) 



Condition (vii): The total electric current Japp is 
known, so the total electric current in material I can be 
written using Eq. (|47]) as 



(B16) 



Ode ~ '^app- 



Equation (B3) then guarantees that the total current in 
material I also equals Japp- 

Condition (viii): There is an arbitrary constant po- 
tential. We set 



Ode 



(B17) 



The eight conditions explicitly given by Eqs. (Bl) 



(B3), (B5), and (B14)-(B17) are general for an isolated 
interface between any two materials I and II, and their 
solution gives the eight unknowns Eq'^^\ Vq^'^^\ V^r*-^'^^^ 

and<"^ 



Oq 



